Quasi-solitons and asymptotic multiscaling in shell models of turbulence 
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A variation principle is suggested to find self-similar solitary solutions (reffered to as solitons) 
of shell model of turbulence. For the Sabra shell model the shape of the solitons is approximated 
by rational trial functions with relative accuracy of O(10 -3 ). It is found how the soliton shape, 
propagation time t n (from a shell # n to shells with n — > oo) and the dynamical exponent zo (which 
governs the time rescaling of the solitons in different shells) depend on parameters of the model. 
For a finite interval of z the author discovered quasi-solitons which approximate with high accuracy 
corresponding self-similar equations for an interval of times from — oo to some time in the vicinity 
of the peak maximum or even after it. The conjecture is that the trajectories in the vicinity of the 
quasi-solitons (with continuous spectra of z) provide an essential contribution to the multiscaling 
statistics of high-order correlation functions, referred to in the paper as an asymptotic multiscaling. 
This contribution may be even more important than that of the trajectories in the vicinity of the 
exact soliton with a fixed value Zo- Moreover there are no solitons in some region of the parameters 
where quasi-solitons provide dominant contribution to the asymptotic multiscaling. 
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I. INTRODUCTION 

The qualitative understanding of many important sta- 
tistical features of developed hydrodynamic turbulence 
(including anomalous scaling) may be formulated within 
the Kolmogorov-Richardson cascade picture of the en- 
ergy transfer from large to small scales. For a dynamical 
modeling of the energy cascade one may use the so-called 
shell models of turbulence which are simplified ver- 
sions of the Navier-Stokes equations. In shell models 
the turbulent velocity field u(k, t) with wave numbers 
k within a spherical shell k n < k < k n+ \ is usually pre- 
sented by one complex function, a "shell velocity" u n (t). 
To preserve scale invariance the shell wave-numbers k n 
are chosen as a geometric progression 



k n = k X r 



(1) 



where A is the "shell spacing" and 1 < n < N. The 
equation of motion reads du n (t)/dt — Q n , where Q n 
is a quadratic form of u m (t) which describes interaction 
of neighboring shells. Clearly, shell models can be ef- 
fectively studied by numerical simulations in which the 
values of the scaling exponents can be determined very 
precisely. Moreover, unlike the Navier-Stokes equations, 
the shell models have tunable parameters (like A) affect- 
ing dynamical features of the of the energy transfer. This 
allows one to emphasize one after another different as- 
pects of the cascade physics and to study them almost 
separately. 

The statistics of u n may be described by the moments 
of the velocity S p (k n ) which are powers of k„ 



S p (k n ) EE (\ Un \ P ) OC h~<* OC A"*- , 



(2) 



in the "inertial range" of scales, ul < n < n>d- Here nr, 
is the largest shell index affected by the energy pumping 
and rid is the smallest shell index affected by the energy 
dissipation. 

In the paper we employ our own shell model called the 
Sabra model IHI. Like in the Navier-Stokes turbulence 



the scaling exponents £ p in the Sabra model exhibit non 
linear dependence on p. Similar anomalies were found 
in the Gledzer - Okhitani - Yamada (GOY) model g],§. 
However the Sabra model has much simpler correlation 
properties, and much better scaling behavior in the iner- 
tial range. The equations of motion for the Sabra model 
read: 

= i(ak n+1 u n+2 u* l+1 + We„"u„ + iu*_i (3) 

a + b + c = . (4) 

Here the star stands for complex conjugation, /„ is a forc- 
ing term which is restricted to the first shells and v is the 
"viscosity" . Equation (|J) guarantees the conservation of 
the "energy" E and "helicity" H 



in the inviscid (y = 0) limit. 



(5) 



In this paper we will consider self-similar solutions of 
the Sabra shell models (||) in a form of solitary peaks 
- solitons. The important role of intense self-similar 
solitons in the statistics of high order structure func- 
tion was discussed in Refs The two-fluid picture 
of turbulent statistics in shell models and corresponding 
"semi-qualitative" theory in the spirit of Lipatovs' semi- 
classical approach |Q was suggested in Refs. pl| , ^2[ : 
self-similar solitons form in and propagate into a ran- 
dom background of small intensity generated by a forcing 
which has Gaussian statistics and ^-correlated in time. 
Accounting in the Gaussian approximation for small fluc- 
tuations around self-similar solitons the authors of pi 12 
reached miltiscaling statistics with a narrow spectrum of 
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z. In the present paper the multiscaling statistics of high 
order correlation functions will be referred to as asymp- 
totic multiscaling. 

Our preliminary direct numerical simulations of the 
Sabra shell model 113 14 shows the asymptotic multiscal- 



ing is a consequence of much reacher dynamics of shell 
models. For example in the 6-interval [— 1 < b < —0.7] 
and at A = 2 we indeed observe extremely intense self- 
similar peaks on a background of small fluctuations. 
Each particular peak has well defined time-rescaling ex- 
ponent zq, however from a peak to peak the value z es- 
sentially varies [I4] . In the region [—0.3 < b < 0] the level 
of intermittency is much smaller, intense solitary events 
vanish, however turbulent statistics remains anomalous 
[14,^6|. Only at intermediate value of b w —0.5 we found 
[13 1 intense self-similar peaks with a narrow spectrum of 
dynamical exponent zq = 0.75 ± 0.02. 

These observations may serve as a starting point in de- 
veloping a realistic statistical theory of asymptotic mul- 
tiscaling which will take into account a wide variety of 
relevant dynamical trajectories of the system not only 
in the vicinity of the well defined solitons. The present 
paper is a first step in this direction and is organized as 
follows. 

The analytic formulation of the problem is presented 
in Sect. 0. For the general reader I describe a self-similar 
f orm of solitary soliton "propagating" through the shells 
( [IB ). I deriv e the "basic self-similar equations" for the 
s olito ns ( II C ) , consider the relevant boundary conditions 
( [ID ) and ana lyze the asymptotic form of soliton tails for 
infinite times (HE). 

Section 111 is devoted to a variation procedure for the 
problem. I suggest in sect. Ill A a simple positive def- 



inite functional T(z) > such that the exact solution 
corresponds to !F{zq) = 0. The a nalytic form of t rial 
functions is discussed in sect. [II B. In section [II C we 
will see in details how the variation procedure works in 
the case of the "canonical" set of parameters A = 2, 
b = —0.5. The characteristic value of J-(z), J-o, is of 
the order of unity. The minimization of T(z) (with re- 
spect of the propagation time and soliton width, with 
proper choice of trial function without fit parameters 
and at experimentally found value zq = 0.75 |l3)) gives 
^min(0.75) f» 0.08. Step-by-step improvement of the ap- 
proximation is reached by a consecutive addition of fit 
parameters which affect a shape of the soliton. With 
10 shape parameters, the value of J- m [ n (z) may be as 
small as 10~ 3 . The resulting "best" values of the shape 
parameters give approximate solutions of the basic equa- 
tion (normalized to unity in their maximum) with local 
accuracy of the order of 10~ 3 . 

Section IV presents results of the minimization on 
trial rational functions (ratios of two polynomials) with 
10 s hape parameters, and their discussion. Firstly, in 
sect. IV A we compare and discuss the z-dependence of 
^minOzTfor b = -0.8, -0.5 and -0.3 (at A = 2). An 



important observation (sect. IV B): there are interval of 
z for which basic equations do not have self-similar soli- 



tary solutions for all times (solitons) but may be solved 
approximately for interval of times from — 00 up to some 
moment in the vicinity of soliton maximum or even after 
it. Configurations of the velocity field in the vicinity of 
these solutions may be called quasi-solitons. My conjec- 
ture is that the quasi-solitons with continuous spectra of 
z may provide even more important contribution to the 
asymptotic multiscaling than the contributions from the 
trajectories in the vicinity of an exact soliton with fixed 
scaling exponent Zq. The concept of quasi-solitons and 
proposed in this paper the dependence of their proper- 
ties on b allows us to reach a qualitative understanding of 
the behavior of intense events for various values of b, ob- 
served in direct numerical simulation of the Sabra model 

In concluding sect. ^ I summarize the results of the 
paper and present my understanding of a way ahead to- 
ward a realistic theory of asymptotic multiscaling for 
shell models of turbulence which also may help in fur- 
ther progress in the description of anomalous scaling in 
the Navier-Stokes turbulence. 



II. BASIC SELF-SIMILAR EQUATIONS OF THE 
SABRA SHELL MODEL 



A. "Physical" range of parameters 



In the inertial interval of scales the Sabra equation of 
motion (||) have formally five parameters: ko, A, a, b 
and c. They enter in the equation in four combinations: 
A, (koa), (ko b) and (fcoc), therefore by rescaling of the 
parameters a, b and c we get ko = 1- Without loss of 
generality we may consider a > 0. A model with nega- 
tive a turns into model with positive (ko a) by replacing 
u n — > — u n . By rescaling of the time-scale t — ► (a t) we get 
a model with a = 1. The arameters a, b and c are related 
by Eq. (Jij). Thus we can express c = — (a + b) — > —(1+6). 
With this choice only two parameters of the Sabra model 
remain independent, b and A. By construction of the 
model A > 1. A typical choice A = 2 will be considered 
in the paper. 

Note that for (a/c) > (which is c > or b < -1 
at a = 1) the model has two positive definite integrals of 
motion which are quadratic in u n : the energy and "helic- 
ity" (H). In this case (as it was discussed in Ref. [^5|) one 
may directly apply the Kraichnan argument for the en- 
stropy and energy integrals of motion in 2D turbulence 
and conclude that in shell models fluxes of energy and 
"helicity" will be oppositely directed: "direct" cascade 
(from small to large shell numbers n) will have integral 
of motion for which large n shells will dominate. There- 
fore for (a/c) > 1 (which is —2 < b < —1 at a = 1) one 
expects direct cascade of "helicity" and inverse cascade of 
energy, like in 2D turbulence. This reasoning predicts di- 
rect cascade of energy for (a/c) < 1 (or b < — 2 at a = 1). 
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In both cases one cannot expect a statistically station- 
ary turbulence with flux equilibrium because one of the 
(positive definite) integrals of motion will accumulate on 
first shells without a mechanism of dissipations. There- 
fore one expects an energy-flux equilibrium with direct 
cascade of energy like in 3D turbulence only for negative 
ratio a/c (or 6 > — 1 at a = 1) when "helicity" integral is 
not positive definite and Kraichnan's arguments are not 
applicable. 

More careful analysis shows that only the region — 1 < 
b < (at a = 1) may pretend to mimic 3D turbulence. 
For b > when (— c) > a = 1 the existence of "helicity" 
integral (even not positive definite) leads to a period-two 
oscillations of the correlation functions which are increas- 
ing with n and "unphysical" from the viewpoint of 3D 
turbulence. One can see this from the exact solution for 
the 3rd-order correlation function 

£3 (&n) = Im(u n _iu„u* +1 ) , (6) 
which in the inertial interval of scales reads || : 

= -lTT; \^e + s(-Y} . (7) 

2k n (a — c) L \as J 

Here e and 8 are fluxes of energy and "helicity" respec- 
tively. For any small ratio S/e and (c/a) < —1 the cor- 
relator Saikn) (and presumably many others) will have 
period-two oscillations which increase with n . Therefore 
in this paper we will consider only the region 

-1< b< 0, 

traditionally keeping A = 2. 

B. Self-similar form of "propagating" solitons 

Self-similarity in our context means that solitons prop- 
agate through shells without changing their form. "Prop- 
agation" means that the time t n at which peak reaches 
its maximum increases with n, in other words, the larger 
n, the later peak reaches this shell. Intuitively this pic- 
ture corresponds to the energy transfer from shells with 
small n to large ones. 

Self-similar propagation of the solitons may be for- 
mally described by the same function f{r n ) of dimen- 
sionless time r„ which is counted from the time of the 
soliton maximum t n and normalized by some character- 
istic time for nth shell T„: 

T n = {t- t n )/T n . (8) 

The time T n has to be rescaled with n as follows: 



, T=l/(k v), 



(9) 



where z is a dynamical exponent and the characteristic 
time T is organized from the characteristic velocity of 



the soliton v and the wave vector of the problem fco . The 
time delays t n — also have to re-scale like T„: 



(10) 



where the positive dimensionless time f is of the order of 
unity. 

Consider next the amplitudes of the solitons. Denote 
as it n! max the maximum of the velocity in the nth shell 
which also re-scales with another exponent y: 



(11) 



In order to relate the exponents z and y we sketch the 
basic equations (^|) having in mind only dimensions and 
A n dependence: 

(du n /dt) oc fco A™m^ . 

Consequently T~ x ~ fcoA™M n!max and therefore 

V + z = 1 . (12) 

Finally we may write a self-similar substitution in the 
form: 



Un (t) = -iv\-y n f{(t~t n )vk \ zn ] 



(13) 



With this choice of a prefactor (— i), the real and positive 
function / will give a positive contribution to the energy 
flux (||) in the inertial interval of scales: 



2fc„ (a - c) Im(u*_ 1 u*ti n+ i) 



(14) 



Note that all these considerations are not specific for 
the Sabra model. They are based on very general fea- 
tures of shell models of turbulence, namely the quadratic 
form of nonlinearity with the amplitudes of interaction 
proportional to A". In Rcfs. |p|,|Tl| similar forms of the 
self-similar substitutions were taken for the Obukhov - 
Novikov and for the GOY models. 



C. Self-similar equation of motion 

Introduce a dimensionless time for the nth shell as fol- 
lows: 



r„ = (t - t n )vk X z 



(15) 



The right-hand side (RHS) of the equation for 
df(T n )/dT n will involve a function / with arguments 
r n ±i and T n ±2- All these times may be uniformally ex- 
pressed in terms of a new dimensionless time 



T EE f/(A Z - 1) 



T n ±s = X ±SZ (T n - T ) +T , S = l, 2 



(16) 

(17) 
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The characteristic time tq is related to the time t n>00 
which is needed for a pulse to propagate from the nth 
shell all the way to infinitely high shells: 



t 



n,oo — / , v m,m-\-l 
m—n 



T \- nz T 



(18) 



Substituting Eq. © in (f§ and using relationship (jl 
one gets the "basic equation" of our problem 



D(r) , "-§1 - C(r) 



(19) 



Here we replaced r„ — > r and introduced a "collision" 
term C(r) according to: 

C(r) = -a A 32 " 2 /* [A 2 (r - r ) + r ] /[A 2z (r - r ) + r ] 

- cA 2 ~ 32 /[A- 2 (r - to) + r ] /[A- 22 (t - r ) + r ] 

+ (a + c) /*[A- 2 (r - r ) + r ] / [A 2 (r - r ) + r ] . (20) 



D. Boundary conditions of the basics equation ( fl9| ) 

The boundary conditions at r = ±oo for a soliton are 
obvious: 



/(±oo)=0 



(21) 



By construction the soliton reaches a maximum at r = 0. 
Therefore 



df(r) 



dT 







(22) 



Introduce a characteristic width of a soliton 1/d accord- 
ing to 



d 2 f(r) 



dr 2 



= -d 2 



(23) 



It is convenient to introduce a new time variable s and a 
function g(s) with the unit width: 



f(r)=g(s), s = rd, 

d 2 g(s) 



ds 2 



= -1 



(24) 
(25) 



With this time variable, the equation of motion ( |19| ) reads 
T>(s, z, s , d) = q(s) - C(s, z, s ) = , (26) 

where 

dg(s) 



ds 



s = rd, 



(27) 



C(s,z,s ) = (a + c) g* ( ' S + s Q ) g [A 2 (s - s ) + s ] 



a A ff *[A 2 (s - So ) + so] g[y z {s - s ) + s ] 



cA 



2-3z 



SO 



+ «o) g 



■ s - Sp 

. X 2z 



+ S ) 



(28) 



Function g(s) should vanish at infinite times: 
lim g(s) = , 

s — ► — oo 

lim g(s) = . 
The boundary conditions at s = read: 



(29) 
(30) 



7(0) = 1 , 



dg(s) 



= 0. 



d 2 g(s) 



s=0 



ds 2 



= 1. (31) 



Note that the problem to find a form of the self-similar 
solitons may be divided into two independent problems: 
for times smaller and larger than so- Indeed, Eq. (^6|) for 
times s < so does not contain functions for times s > so 
and vice versa. At the boundary between these regions 
Eq. (p6|) reduces to 



d r l 

dso 



g(so) 



dpo , 



(32) 



po = a(X Az - z - 1) - c(l - \ z - iz ) 



In our discussion a > 0, c<0, A>1 and z > 2/3. 
Therefore po > and dg(so)/dso < 0. It means that the 
time so is larger than the time s = at which g(s) has a 
maximum, i.e. s > 0. As we noted [and see Eq. (jig)], 
so is the time which is needed for a pulse to propagate 
from the nth shell to infinity high shell. Therefore the re- 
lation so > agrees with our understanding of the direct 
cascade. 

So, we will divide the time interval [— oo < s < oo] into 
two subintervals: [— oo < s < so] and [so < s < oo]. For 
[— oo < s < So] the maxima of the solitons in all shells 
are in the inertial interval of scales and very high shells 
are not yet activated. We will refer to this interval as 
an inertial interval of times. In the second time interval, 
for s > so high-shell solitons already reached the dissi- 
pative interval of scales. We will refer to this interval 
as a dissipative interval of times. Generally speaking, in 
this interval one has to account for the viscous term in 
the equation of motion. In this paper we will restrict 
ourselves to the inertial interval of times. 

It was shown in Ref. || that equations similar to ( JTs| ) 
with similar boundary conditions can be considered as 
nonlinear eigenvalue problems. They have trivial solu- 
tions /(s) = 0, but they may have nonzero solutions for 
particular the values of z = zo and so- Below we will find 
nontrivial solutions of our Eq. (|2^) that satisfies condi- 
tions (09) , (plh and for which zo lies in the physical region 



< z < 1. 
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E. Qualitative analysis of self-similar solutions at 

s — > ±oc 

In this subsection we will analyze time dependent solu- 
tions of (||) which are more general than a Kolmogorov-41 
(K41) and have a form of solitary pulses - solitons. We 
will show that the solitons have long (power-like) tails 
g(s) oc s x± at s — > ±00. One way to find the asymp- 
totic solution of Eq. ( |26| ) is to balance the exponents in 
its left-hand side (LHS) and RHS. For g(s) oc s~ x this 
gives immediately x = 1. Thus: g(s) — Di/s. Equating 
prefactors in the LHS and RHS of the Eq. (|l9|) we get: 



D 1 = 



-X 2 d 



{a-c\ 2 ){\ 2 - 1) 



(33) 



The coefficient D\ appeared to be real and, in the actual 
range of the parameters (A > 1 a > 0, c < 0), nega- 
tive. For the positive function g(s) this asymptotic form 
describes the front part of the pulse (at negative s). So: 



g(s) 



-X 2 d 



s(a-cA 2 )(A 2 - 1) 



-00 



(34) 



Another approach is to assume that in the equation 
[g(s) oc s~ x ] the exponent x < 1. Then, the LHS will 
behave as l/s 1+x while the RHS will be proportional to 
l/s 2x . In the limit s — ■> 00 and at x < 1 one may neglect 
the LHS of the Eq. (|2^). Then the exponent x cannot be 
found by power-counting. Instead one requires that the 
prefactor in the RHS [i.e. in the C(s) term] must vanish. 
This gives the following equation for x: 



= aA 3 ^ 1 -*)- 2 + cA 2 - 3 ^ 1 "") - (a + c) 



(35) 



Denoting 

A = , (36) 

we have, instead of Eq. (|35|), the square equation for A: 
ah 2 - (a + c)A + c = , 

with the roots 

Ai = 1 , A 2 = c/a . 

In the chosen region of parameters (c and a have different 
signs) A2 < 0, which contradicts the assumption of real 
x. Therefore the only root Ax = 1 is relevant. According 
to definition (|36]) this gives 



xi = 7 



1 



2 

3z 



(37) 



Here we used notation 7 for the exponent of the long pos- 
itive tail of the pulse introduced in Q for the Novikopv- 
Obukhov shell model. As we see, the relation ([3?]) is 
model independent. Actually this equation is a conse- 
quence of the conservation of energy (reflected in the 



constraint a + b + c = 0) and the fact that the shell 
models account only for interaction of three consequent 
shells. Clearly, the relevant region of 7 is < 7 < 1. 
This corresponds to 



-<z<l 



(38) 



We conclude that for s 



g(s) = D, ( / S i , (39) 

with a free factor D y which has to be determined by 
matching the asymptotics (39j) with a solution in the re- 
gion s ~ so. It is known p] that the self-similar core of 
the function g(s) gives rise to a linear n-dependence of 
the scaling exponent of the nth order structure func- 
tion at large n: 



( n = z + (1 - z)n , at n > n CI = I/7 



(40) 



For 2=| Eq. (|39| ) gives the K41 slope | which is 
the largest possible one. Value z = 1 corresponds to the 
largest possible intermittency (zero slope) . Therefore we 
may consider Eq. ([38]) as the "physical region" of the 
dynamical exponent z. 



III. VARIATION PROCEDURE FOR SOLITONS 
AND QUASI-SOLITONS 

A. Suggested functionals 

Consider the following functional acting on the func- 
tion g(s) 



F{g(s)\z,s ,d} = 




ds[D(s, z, so, d)}' 2 



(41) 



which also depend on the parameters of the problem 
A, b, z, so and d. By construction, the functional 
T{g(s)\z, so, d} is non-negative and equal zero if g(s) is 
a solution of the problem (p6|): T> = 0. Clearly, the func- 
tional ( [fl| ) is not unique. For example, we may use more 
general, "weighted" functional, 



F{g{s)\z, s , d} 




dsW(s)[D(s,z,s ,d)} 2 , (42) 



with some positive "weight" function W(s) > 0. This 
functional is also positive definite and also equal zero if 
g(s) is a solution of the problem (^6|). 

One can easily find many other functionals giv- 
ing approximate solutions of the problem. The func- 
tional J 7 {g(s)\z, so, d} has an advantage of simplicity. Its 
minimization (with a proper choice of the trial func- 
tion) leads to a solution of the problem with hight ac- 
curacy. For example, with trial functions discussed in 
Section [II B the relative accuracy (with respect to the 
value of a soliton maximum) is O(10~ 3 ). Therefore we 
will restrict the present discussion to the simple func- 
tional J-~{g(s)\z, sq, d}. 
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B. Suggested form of trial functions 

For simplicity, in this paper we will seek only real 
(without nontrivial phases) solitons. Complex solitons 
will be discussed (within the same scheme) elsewhere. 

My suggestion is to use different trial functions for neg- 
ative and positive times, both satisfying boundary con- 
ditions ( J3l| ) at s = 0. Denote as g m+ (s) the trial func- 
tions for positive times with m shape parameters. Let 
the function g m+ (s) satisfy condition ( [32] ) at s = sq. 
This condition is a constraint on the time derivative of 
l/<7 m +(s). Having also in mind that g m +(s) is defined 
on the interval [0 < s < Sq] with s ~ 1 it is convenient 
to choose l/g m +(s) as a polynomial in s. The first (free) 
term of expansion is 1 because g m +(0) = 1. The sec- 
ond term (oc s) vanishes due to dg(s)/ds — at s = 0. 
The next term must be s 2 /2 due to the restriction on 
the second derivative (^TJ) at s = 0. In order to satisfy 
condition ( [32] ) at s = sq we have to account at least for a 
cubic term. Therefore the simplest function of this type, 
<7o+(s) without fit parameters takes the form 



90+ 00 = 



1 



3s 2 



d 



so 



(43) 



If needed we can add fit parameters pi, p 2 , etc. account- 
ing for terms cx s 4 , oc s 5 , etc. Accordingly, a trial func- 
tion with m shape parameters takes the form: 



9m+(s) = 



2 3 m 

1 + TT + 7T9 ( - s - ^Pj) 

3=1 



+£77 



3sg V d 

«J+3 



In most cases it would be enough to use g^ + (s) with 5 
shape parameters. 

Denote as <? m _(s) the trial functions for negative times 
having m shape parameters. These functions have to ap- 
proximate the basic equation on the infinite time interval 
from — oo to 0. Therefore their analytic form is a much 
more delicate issue. Analyzing the form of the peaks ob- 
served in direct numerical simulations |p^| . I have found 
a function gi-(s) with one fit parameter q, 



[1 - s/y/qfa-l)] 



(44) 



which allows us to reach accuracy of solution of Eq. (|2£ 
of 0(1CP 2 ). This accuracy would be enough for future 
comparison of the "theoretical" shape of a pulse with that 
found in direct numerical simulation of the Sabra model. 
However this function does not have the correct asymp- 
totic behavior ( p4[ ) for s — > — oo and is inconvenient for 
successive improvements of the approximation. 

Instead of <7i-(s) f° r negative time I use in the reg- 
ular minimization procedure a ratio of two polynomials 



of nth and n + 1 orders. This ratio has 2n + 2 free pa- 
rameters. After accounting for three conditions (|3l]), just 
(2n — 1) parameters remain free. Simple function of this 
type with n = 1 , which also agrees with asymptotic (|^|) , 
has no free parameters: 



§o-00 



D lS 



(45) 



We will discuss even more simple rational function 

50 -(.s) = [l + i S 2 ]- 1 . (46) 

Actually a better approximation may be achieved by a 
function similar to ([45|) with one free parameter qi : 



9i-( s ) 



1 + q lS 



l + q lS 



1*2 • 



(47) 



In actual calculations it would be sufficient for our pur- 
poses to account for 5 shape parameters in the function: 



35- (s) 



l + qis + q 2 s 2 



g 4 s 



1 + q lS + (q 2 + i)s 2 + q 3 s 3 + q 5 s 4 



(48) 



The function # 5 _ (s) reduces to 33- (s) by choosing (74 
q 5 = and to gi-(s) at q 2 = q 3 = <? 4 = <7s = 0. 



C. Test of the variation procedure 

In this subsection we will see in detail how the vari- 
ation procedure suggested above works for the "canon- 
ical" set of parameters A = 2, b = —0.5. Denote as 
Fm\n,m{z) the result of a minimization of the functional 
J 7 {g(s)\z, sq, d} acting on the functions g m ±(s) at given 
z. The parameters in the minimizations are: the propa- 
gation time so, width parameter d, m shape parameters 
for negative times (31, •• • q m ) and m shape parameters 
for positive times (p±, ■ ■ ■p m ): 



{z)= ruin T{g(s)\z, s ,d} 

Qj,Pj\ s o,d, 



(49) 



Table 1 displays the found values of J r m in,m(-2 ; ) for z 
0.75. 0.85, 0.95 and m = 0, 1, 3, 5. 

Consider first the case z = 0.75, which corresponds 
to the dynamical exponent ofintense self-similar peaks 
observed in our direct numerical simulation of the 
Sabra shell model Jl^]. At m — 5 we have reached 
^"min,5(0.75) w 2.3 • 10~ 4 . This value is more than 200 
times smaller than characteristic value of the functional 
before minimization, J-q w 0.5. This allows us to hope 
that the minimization using the functions (?5±(s) with 
total number 2 m = 10 of the shape parameters will be 
sufficient for most applications. 

Note that for some purposes we may use less then 10 
shape parameters. For example "optimal" values of sq 
and d shown in the Table 1 begin to converge for m > 3. 
Therefore for a reasonable good estimate of Sq and d we 



G 



may use the trial functions <73±(s) having 6 shape pa- 
rameters. Moreover, the same level of accuracy as with 
772 = 3 may be achieved just with two shape parameters 
(q, pi) if we replace g%- (s) by gi-(s) and, (which is less 
important) 53+ (s) — > (Ji+(s). Minimization with <?i_(s) 
and (s) gives: 

min F{g{s)\z, s , d} ps 0.0197, at z = 0.75, (50) 

s ,<2,i},pi 

s ps 0.54 , d ps 1.57 , g ps 8.59 , pi ps -0.44 . 



Another example in which we may use less then 10 
shape parameters for a reasonable good description is the 
shape of solitons, g(s). Fig. |l| displays functions g(s) and 
q(s), Eq. ( p7j ) for 777 = 3 (dashed lines) and for m = 5 
(solid lines). We see that for comparison with experi- 
ments we can use simpler form with m = 3. In Fig. |l| 
are also plotted by dashed lines functions g(s) and q(s) 
resulting from minimization ( p0[ ) with 2 shape parame- 
ters. They are indistinguishable from the corresponding 
functions having 6 shape parameters. 

Table 1. Optimal values of (sq, d) for z = 0.75 and 
values of J- m in,m for z = 0.75, 0.85 and 0.95. Number of 
the shape parameters m = 0, 1, 3 and 5. 
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FIG. 2. Local deviation X> m i n ,5(s, 0.75), Eq. ( pi| ) vs s for 
A = 2, b — —0.5. Upper panel: region —10 < s < so. Lower 
panel: blow up of the region —so < s < sq. 
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FIG. 1. Trial functions g(s) and q(s) = d[d(;(s)/ds] for 
A = 2, 6 = —0.5 and 2 = 0.75. Solid lines correspond to 
the full minimization (m = 5), dashed - to m = 3. The same 
dashed lines show results of the minimization (^FJ) with only 
two shape parameters. 



Note that the minimal values of the functional 
Fmin,m{ z ) characterize the "global" accuracy of the ap- 
proximation in the whole interval of s, [00 , sq]. It would 
be elucidative to discuss a "local deviation" which may be 
described by the function V(s, z, So, d) = q(s)—C(s, z, So). 

According to Eq. (p6|), the local deviation "D(s, z, s , d) 
has to vanish for all s. Denote as P m i n ,m (s, z) the de- 
viation V{s, z, so,d) at optimal values of so,d and 2m 
shape parameters. The deviation 2? m i n) 5(s, 0.75) is shown 
in Fig. H Excluding the small region |s| < 0.1 the de- 
viation £> m in,5(s, 0.75) < 10 -3 , i.e. in more than 500 
times smaller than the maximal value of q(s) shown in 
Fig. |l|. This serves for us as a strong support of the 
conjecture that by adding more and more fit parameters 
one can reach smaller and smaller values of ^"min, i.e. 
lim„ woo J- m in.m = and for z — 0.75 one can find a true 
solution of the problem. The currently available personal 
computers allow to find solutions with .Fmin.io ~ 10~ 3 
and V < 2 • 10~ 3 during one-two hours of calculations 
using the standard package of Wolfram' Mathematica. 

Consider now different a value of z — 0.95. As fol- 
lows from the Table 1, there are no jumps in improving 
the approximation for m > 2. For z — 0.95 the ratio 
3~ m i n ,i(z) / Tmin^z) ~ 2 while for z = 0.75 the same ra- 
tio J 7 mm.i(z)/J-mhi,b(z) — 28. It is very reasonable to 
expect that ^ r m in,m(0.95) cannot be much smaller that 
0.06 even for very large m. One may conclude that for 
z = 0.95 there is no "global" solution in the interval 
[— oo < s < sq]. 
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FIG. 3. Local deviation D m m,5(s, 0.95) vs s for A = 2, 
b — —0.5 and z = 0.95 with optimal set of 10 shape pa- 
rameters. Upper panel: region — 10 < s < so. Lower panel: 
blow up of the region —so < s < so. 

Nevertheless, the local deviation T> m i n ^(s, 0.95) shown 
in Fig. |^ is quite small, say in the region s < 0. Moreover, 
in the region s < s» « — 2 the characteristic value of the 
deviation V min ^(s, 0.95) is more or less the same as for 
z = 0.75 where we have approximate solution (soliton) in 
the whole region s < sq. Consequences of this fact will 
be discussed later. 



IV. SOLITONS, QUASI-SOLITONS AND 
ASYMPTOTIC MULTISCALING 



A. z-valleys of the functional J r mm,s(z) 



In this section we discuss the ^-dependence of the func- 
tional ^min.s (z) and of the optimal values of sq and d. 
These functions for b = —0.5, —0.8, —0.3 are displayed 
in Figs. ||,||, ^ respectively. 

For b = -0.5 (Fig. |) the z-dependence of the func- 
tional has a minimum J- m i n ,5(z) ~ 2 • 10 -3 around 
z w 0.734. This minimum is quite flat. For example, 
^ 7 min,5(^) < J~* = 2.5 ■ 10 -3 for z within the interval 
[zmin = 0.71, z max = 0.76]. Note that for b = —0.8, 
(Fig. ||) the same (arbitrary) level exceeds T m - in ^(z) 
for z in the wider interval [0.71,0.84]. In the same time 
for b = -0.3 (Fig. |) there arc no values of z for which 
J~ m i n £(z) < T*. A natural interpretation of these facts 
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dashed lines shows K41 value of z = 2/3. Horizontal dashed 
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Points for 10 3 ^ min b with 0.72 < z < 0.76 are below this level. 
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FIG. 5. The z dependence of 10 3 < F, n i» l 5, 10 2 ^" m in,5 (Upper 
panel) and so, d (Lower panel) for A = 2, b — —0.8. Vertical 
dashed line shows K41 value of z = 2/3. Horizontal dashed 
line in upper panel corresponds to chosen value 10 3 JF, = 2.5. 
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vertical dashed line shows K41 value of z = 2/3. Horizon- 
tal dashed line in upper panel corresponds to chosen value 
10 3 JT» = 2.5. All points for 10 3 J-" m in,5 exceed this level. 

is that at b = —0.8 one can meet in the velocity real- 
ization (of the Sabra model) u n (t) intense solitons with 
values of z in a quite wide interval (for example in the 
interval [0.71, 0.84] mentioned above), at b = —0.5 "al- 
lowed" interval of z is more narrow (say, [0.71, 0.76]) 
and for b = —0.3 one hardly can meet intense solitons 
at all. As we will discuss below this statement is in a 
quantitative agreement with the preliminary results of 
the numerics ]i"4j . 

Clearly, we are not talking about particular values of 
z m m and z max , for example because the boundary level 
= 2.5 • 10~ 3 was chosen arbitrarily. Moreover, the 
objects z m j n , z max do not have explicit sense. The actual 
conjecture is that the functional J- m i n , m {^, b, z) is corre- 
lated with a probability to meet a soliton or quasi-soliton 
with given dynamical exponent z in the realization (for 
given A and b): the smaller value of J- m i n ,m(z) (at large 
enough m), the larger this probability. From this point of 
view the two following scenarios are statistically almost 
equivalent. The first one may be called a multi- soliton 



scenario: there is a discrete spectrum of solitons with 
some close set of exponents z 0j i, £0,2 m the interval [z m j n , 
ZmaxL, O ne can imagine that this is the case by looking at 
Fig. [5] (upper panel) where the function F m - m ^(2, —0.8, z) 
has two minima at zo,i ~ 0.77 and 20,2 ~ 0.81. The 
second one will be referred to as a quasi-soliton scene- 
rio with continuous z-spectrum of quasi-solitons in some 
(wide) interval of z. 



B. Local deviations and quasi-solitons 



The analysis of the local deviations T> m i lhm (s, z) pre- 
sented in this section supports a quasi-soliton scenario of 
asymptotic multiscaling. Let us chose b = —0.8 for which 
we found in Fig. (upper panel) the wide deep z- valley 
of the functional .Fmin.s (z). Compare the s-dependence of 
Anin,m(s, z) for the value of z = 0.81, which corresponds 
to the right local minimum of the functional and for two 
slightly larger values of z = 0.84 and 0.88. The local de- 
viations X> m in,m(s, z ) are plotted in Fig. Q for z = 0.81 by 
solid lines, for z = 0.84 by dashed lines and for z = 0.88 
by dash-dotted lines. The upper panel shows s-interval of 
a left tail the soliton up to its maximum: [—10 < s < 0]. 
The lower panel shows s-interval around soliton maxi- 
mum: [—1 < s < sq]. 
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It was mentioned that the value zo,2 = 0.81 corre- 
sponds to the right local minimum of the functional. 
The local deviation \T> m i n ^(s, 0.81) | is about 2 • 10~ 3 
in the s-region [— l,So], see solid line in Fig. [t], lower 
panel). For s < — 2 (the upper panel) the local devia- 
tion is almost 10 times smaller, about 2.5 • 10~ 4 . The 
value z = 0.84 does not correspond to z-minima of the 
functional: J c mi „, 5 (0.84) « 2 J" mi „, 5 (0.81). The larger 
value of ^ r min,5(0.84) is mostly determined by the re- 
gion of positive s where \V m i ny5 (s, 0.84) exceeds 0.01, 
which essentially larger than the local deviations in this 
region at z = 0.81, [compare dashed (z = 0.84) and solid 
(z = 0.81) lines in the lower panel of Fig. [7J. Unex- 
pectedly, in the region s < — 1 (the upper panel) the 
local deviation at z = 0.84 even smaller than that at 
z = zq ^ = 0.81. In other words, in this s-region the 
reached approximation is better for value z = 0.84 than 
that for z — 0.81. Dash-dotted lines in Fig. ^| show that 
the local deviation for a z = 0.88 essentially larger than 
for z = 0.84 and 0.81. 

Let us show that these qualitative results are indepen- 
dent of the particular form of the functional (|4l]) which 
we chose. For that consider the weighed functional ji^ ) 
with the weight fuction: 



W(s) = [(s - SQ ) 2 + If : 



(51) 



which is of the order of unity in the vicinity of the soliton 
maximum and increases toward the front tail as \s\ 2y . At 
y = one recovers the uniform weight W{s) = 1. For 



positive y the weight function emphasize the left tail, 
i.e. the region s < 1. The value of y has to be smaller 
than ~, otherwise the s- integral in (^) diverges. It is 
reasonable to choose an intermediate value of y = 1/4 
and to repeate the minimization procedure. The result- 
ing changes in the values of sq, d and shape parameters 
were minor. For example (for A = 2, b = —0.8, z — 0.84) 
s = 0.7506 -> 0.7478, d = 1.2645 -> 1.2618. Fig. | com- 
pares the local deviation V for unweighed (solid lines) 
and weighed (dashed lines) functionals. Again, the dif- 
ference between these two approaches is minor. 

A natural interpretation of all these facts as is follows: 
for some discrete values of z = zqj (say for zo,i ~ 0.734 
at b = -0.5 and z ,i ~ 0.78, z 0:2 ~ 0.81 at b = -0.8) 
the variation procedure allows to find better and bet- 
ter approximation to the solution in the global sense: 



in the limit m — > oo. Clearly in these 



cases solution g{s) exists also locally: the maximum of 
the local deviation T> m i n rn (s, z j) also goes to zero in 
this limit. In means that there are self-similar solitons 
in the shell model with scaling exponents z 0j -. For other 
values of z there is no global solution. It means that 
for to — > oo limit of J- m in.m is finite, but may be small 
enough. Nevertheless for a wide region of s < s* (in 
the example z = 0.75, b = —0.5 the value of s* w —2) 
the local deviation P m in,m(s, z) is very small and equa- 



tion (|26|) may be "almost solved" . We may say that 
the variation procedure finds "quasi-solutions" , which 
approximate the equations of motion with very good ac- 
curacy starting from the "time of appearance of a soliton 
(s — > — oo) until some time s*. The value of s* may 
be close to the time of the soliton maximum. Remem- 
ber, that in the real shell systemintense events form in a 
random background of small fluctuations where one may 
meet configurations corresponding to initial conditions of 
these "quasi-solutions" . In that case the system begins to 
evolve along the quasi-solutions up to some time about 
s < s* where the quasi-solution stops to approximate 
well the self-similar equation of motion. These peaces 
of trajectory are referred to as quasi- solitons. After the 
time s* trajectory has to deviate from the quasi-solitons, 
presumably exponentially fast. Quasi-solitons look like 
the front part of solitons and may include their maxima. 

The main message is that the quasi-solitons have a con- 
tinuous spectrum of z which vary a lot around zq . There- 
fore their contribution to the asymptotic multiscaing may 
be even more important than the contributions of trajec- 
tories in the vicinity of the soliton with fixed scaling ex- 
ponent zq- Moreover in some region of parameters [in 
Sabra model at A = 2 for b > b CI « -(0.5 - 0.4)] there 
are no solitons and the quasi-solitons provide dominant 
contribution to the asymptotic multiscaling. 

This picture is consistent with the preliminary direct 



numerical simulation of the Sabra shell model [|14 13 
where it was observed: i) at b = —0.8 self-similarintcnsc 
events with different rescaling exponents (see wide deep 
minimum of ^ r m in,m(^) in Fig- ^); ii) at 6 = —0.5 self- 
similar events with a very narrow region of z (see not so 
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deep minimum of J-~ m in,m{z) in Fig. ^); and iii) no solitons 
at b = —0.3 (there is no deep minimum in Fig. ||). 



V. CONCLUSION 

• I have suggested a variation procedure (basic func- 
tional and trial functions) for an approximate solution of 
equations for the self-similar solitary peaks (solitons) in 
shell models of turbulence. The rational functions with 
10 shape parameters approximate solitons in the Sabra 
shell model with relative accuracy 0(10 -3 ). 

• For the standard set of parameters (A = 2, b = —0.5) 
the dynamical exponent zq rs 0.734 found in the paper 
agrees within the error bars with the experimental value 
of z w 0.75 ±0.02 fl|. 

• The variation procedure allows to find trajectories of 
the system which are very close to the self-similar soli- 
tons during an interval of time from — oo up to some 
time in the vicinity of the soliton maximum or even after 
it. These trajectories, called quasi-solitons have contin- 
uous spectrum of the dynamical scaling exponents and 
may provide dominant contribution to asymptotic mul- 
tiscaling. The discovered features of quasi-solitons for 
b = —0.8 , —0.5, and —0.3 at A = 2 allows me to ra- 
tionalize the preliminary numerical observations [ p~3|Jl4) 
of asymptotic multiscaling for various values of b in the 
Sabra model. 

• This paper may be considered as a first step toward 
a realistic statistical theory of asymptotic multiscaling in 
shell models of turbulence which will account for a wide 
variety of trajectories of the system in the vicinity of 
quasi-solitons. In particular, one can apply the variation 
procedure to find possible complex solitons and quasi- 
solitons with non-trivially rotating phases. Such objects 
may be important in quantitative description of high- 
level intermittency when the parameter b approaches the 
critical value —1. Analytical form of self-similar trajec- 
tories of the system found in this paper may help one 
to analyze stability of nearby trajectories of the system. 
This is useful for describing how the system approaches 
and later escapes a vicinity of quasi-solitons in order to 
describe their role in statistics ofintense but not solitary 
events. 

I hope that the analysis of dynamics ofintense events 



in the Sabra model, presented in the paper will help in 
futher understanding of asymptotic multiscaling along 
similar lines. These will be also useful for further progress 
in the problem of anomalous scaling in Navier-Stokes tur- 
bulence at least on a qualitative and may be on semi- 
quantitative level. 
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